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' Abstract 

^ |. A Fokker-Planck equation on fractal curves is obtained, starting from 

I ' Chapmann-Kolmogorov equation on fractal curves. This is done using the 

' recently developed calculus on fractals, which allows one to write differen- 

ce , tial equations on fractal curves. As an important special case, the diffusion 

and drift coefficients are obtained, for a suitable transition probability to 
get the diffusion equation on fractal curves. This equation is of first order 
in time, and, in space variable it involves derivatives of order a, a being 
the dimension of the curve. The solution of this equation with localized 
^ , initial condition shows deviation from ordinary diffusion behaviour due to 

underlying fractal space in which diffusion is taking place. An exact solu- 
tion of this equation manifests a subdiffusive behaviour. The dimension 
of the fractal path can be estimated from the distribution function. 

Fractal curves and paths are encountered frequently in Physics [Ij [2]. Sev- 
' eral geometries like that of polymer chains, percolating clusters, Brownian and 

[ Fractional Brownian trajectories and many more have been identified as fractals, 

more precisely as fractal curves [T] . Transport on these structures reveal remark- 
able properties [31l31[51[Sl[71IHl[n]- In particular, anomalous diffusion on fractals 
■ is a topic of immense current interest [TUl[TT 1 [n i [n i [TIl[T51[TB l [r7 l [T51[Tg]. 

I There have been several investigations [20l [211 HH [23] , both analytical as well 

' as numerical which shed light on the various facets of the problem. These in- 

vestigations use the Fractional operators . The Fractional Derivatives are often 
used to explore the charcateristic features of anomalous diffusion by setting up 
fractional kinetic equations. Fractional derivatives are non local operators and 
hence not always suitable to handle local scaling behaviour. Analysis on frac- 
tals is another remarkable development which has been extensively used for the 
treatment of diffusion, heat conduction, waves etc on fractals. But, in this ap- 
proach the operators are constructed using the self similarity of fractal sets, and 
these are restricted to such sets which are post critically finite. However, as is 
well known, an ordinary diffusion equation is inadequate to describe anomalous 
transport, and an exact equation based on an appropriately developed calculus 
is still desired. 
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Ordinary calculus does not equip us to handle problems such as anomalous 
diffusion, dynamics on fractals, fields of fractally distributed sources etc by 
setting up and solving ordinary differential equations. Several authors have 
recognized the need to use fractional derivatives and integrals to explore the 
characteristic features of fractal walks, anomalous diffusion, transport etc. by 
setting up fractional kinetic equations, master equations and so on. 

Further measure-theoretical approaches are used which include defining deriva- 
tive as inverse of the integral with respect to a measure and defining other 
operators using derivatives. Even though measure theoretical approach is ele- 
gant, Riemann integration like procedures have their own place. They are more 
transparent, constructive and advantageous from algorithmic point of veiw. 

Fractal curves lack the smoothness properties required by ordinary calculus. 
Recently, [24] a new calculus on fractal curves, such as the von Koch curve, was 
formulated. In this calculus, a Riemann-like integral along a fractal curve F is 
defined. This integral is called the F"-integral, where a is the dimension of the 
curve. A derivative along the fractal curve, called F"-derivative, is also defined. 
These operators are different from the fractional operators used in [TUl El [T^] , 
in the sense that they are local and a newly defined measure like quantity called 
the mass function (which is algorithmic in nature) is used to define these. The 
order of these operators is exactly equal to the dimension of underlying space 
, and thus they reduce to ordinary integral and derivative operators when the 
space is R and not a fractal. 

Several aspects of this calculus retain much of the simplicity of ordinary 
calculus. In fact a conjugacy between this calculus and ordinary calculus on the 
real line is established. This new calculus which is simple, direct and algorithmic 
can be applied to various physical processes. 

Here, using the framework of this calculus, we develope a first-principles 
approach to the Fokker-Planck equation on fractal curves. A particular choice 
of transition probability then leads to a new form of diffusion equation. This 
equation is of the first order in time and involves application of the i^"-derivative 
with respect to spatial variable twice. An exact solution of this equation shows 
a subdiffusive behaviour. 

We begin by fixing our notation as in [21]. We consider a (fractal) curve 
F C R" which is continuously paramatrizable i.e there exists a function w : 
[ao,&o] F C R" which is continuous. We also assume w to be invcrtible. 
A subdivision P[a,b] of interval [a, 6], a < 6, is a finite set of points {a = uq < 
ui, . . . < Un — b}. For gq < a < b < bo and appropriate a to be chosen, let 



where | • | denotes the Euclidean norm on R". and \P\ = max{uj+i — Ui\i = 
0, . . . , n — 1} A new dimension, the 7-dimension of F , which will be denoted by 




(1) 



dim-y(_F), is given by 



dim.^(F) 



inf{a : 7"(i^, a, 6) = 0} = sup{a : 7"(i^, a, b) 
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Hereafter, a will be assuraed to be equal to dim^(i^) (thus, a > 1). 
The rise (staircase) function 5"^ : [ao,6o] — R of order a for a set F, is 
defined as 

where < pq < bo is arbitrary but fixed, and u e [ao,6o]. It is a monotonic 
function. We denote a point on the fractal curve F by 9 — w(u), and define 

J(0) = 5^(w-i(0)), eeF 

Hereafter we consider only those curves for which Sp is finite and invertible 
on [a,b]. The F"-derivative of a bounded function f : F ~^ li a.t 9 £ F is 
defined as 

(Dlfm) = F- lim ■^(^'j ~ ^^^\ (3) 

where the F- lim denotes limit along points of F, if the limit exists. 

Let C(a, &) denote the segment {w(it) : u G [a, 6]} of F. A Riemann- 
like integral on F, called F"-integral, is also defined [24|- It is denoted by 
/c(a 6) /(^)'^F^- The above mentioned F"-derivative and F"- integral are related 
to each other through the Fundamental Theorems of Calculus as "inverses" of 
each other [53] . 

The notion of conjugacy of calculus on fractals and ordinary calculus on 
the real line is very useful. Let denote the map (conjugacy) from the class 
of bounded functions on F to the class of bounded functions on the interval 
[5^(ao),5'^(6o)]defined by 0[/](S'^(u)) = /(w(u)) Then it follows gl] 

g{x)dx^ / .f{9)d%9 (4) 

S°(a) Jc(a.b) 

where g = '/'[/]• Moreover (j) relates derivatives Dp with the ordinary deriva- 
tive D, thus {Df.f)yv{u) = {Dg){u). 
The Taylor series is given by 

h{9) = y i:li^l—Wl{D'^rh{9') (5) 

provided the bounded function h is F"- differentiable any number of times on 
C'{a,b). It is also possible to write a Taylor series with remainder. 

Let J^^^ V{9, t)dp9 be the probability that a particle constrained to move 
on F is found in the segment C{a,b), or in other words, V{9,t) denotes the 
'fractal' probability density that the particle is found at 9 at time t. Let the 
probability density for transition from a point 9' at time t to at time t + r, be 
denoted by P{9,t + T\9' , t). A formalism to analyse similar situation in ordinary 
space is developed in ,25 , we intend to modify the same for the case of fractal 
curves. 
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The Chapmann-Kolmogorov equation on fractal curve F can be written in 
the form 

Vi9,t + T)= [ P{9,t + T\e',t)V{e',t)d'^9' (6) 

JC{a,b) 

where 9,9' € F. Let A = A{9,9') = J{9) - J{9'). The intergrand in equation 
® is: 

P{e,t + T\9',t)V{9',t) = P{J-^iJ{9)^A + A), 

t + T\J-\j{9) - A),t) 
xV{J-\j{9)-A),t) 

The Taylor expansion of this integrand in ^ then leads to 

V{9,t + T) - V{9,t) 

/°° (—^^^ 
E S^A"P^|,)"{P( J-i(J(0) + A), 
n=l 

t + T\9,t)V{9,t)}d1^9' (7) 

While there are other ways to write the Fokker-Planck equation we turn to 
the use of conjugacy which reduces the problem to the ordinary case and the 
meaning of moments will be transparent. We use the following explicit notation 
for the conjugacy map cf) : 

{cj,gH){Ji9),9')=H{9,9') 

{<l>e'H){9,J{9')) = H{9,9') 

while applying to a function H{9, 9') of two arguments 9 and 9' . 

Let y = J{9),y' = J{9') and assume S'^(ao) <y,y'< Sf{bo). Further let 
us denote 

A{9,y')^q^g,A{9,9') 
A'{y,y')^^gA^<f>eo(j,e'A{9,9') 

Then, 

A{9, y') = J{9) - y' and A'(y, y') ^ y - y' 

where S^ao) < y,y' <Sf (bo). 

Now we define V — (j)e{V) and P' = (j)0> o (j)gP. Using the conjugacy of 
integrals from equation equation Q becomes: 

r{y,t + r)^V'{y,t) = E / 



i^r{P'iy + A',t + r\y,t) V'{y,t)}dy' (8) 
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Integrating over A', we see that dy' = —dA' Hence, 

V'{y,t + r)-V'iy,t) = - ^ L_L / (A^ 



n— 1 



y-S°(ao) 



{^r{P'{y + A',t + T\y,t) V'{y,t)}dA' (9) 

dy 

The transitional moments are given by 



Mn{y,t,T)^(f>e'Mni9\t,T) (10) 

ry-Sp{ao) 

{ATP'{y + A',t + T\y,t)dA' (11) 



y-S°(fco) 



Hence 



M^{9',t,T)= {j{e)^ j{e')YP{e,t + T\e\t)d%e (12) 

"'C(ao,bo) 

Substituting PU)) in equation (|n]),and applying conjugacy we get 

v{e,t + r) - v{e,t) ^ f^i-D-^r \ ^^"^^f^^ v{0,t) 

Now we assume that the moments Af„ and M„ can be expanded in a Taylor 
series. 



= iAr>{y,t)T + OiT') (13) 

It: 

and 



Af„(M,r) ^^(n)(g^^)^^o(^2) (14) 

The term of order r" vanishes because for t = the transition probability 
is P' {y,t\y' ,t) = S{y,y') which leads to vanishing moments. By taking into 
account only the linear terms in r we have the Kramers-Moyal expansion 

rj OO o 

g-V'{y,t) = Y^i- r{(Ar)(y,t)V'{y,t)} (15) 

n=l ^ 

i.e. 

^vi0,t) = 5](-^?^)"{A(")(^^,^)v^(^?,^)} (16) 

n— 1 

From the above we obtain the Fokker -Planck equation if the expansion in 
equation ([TS]) stops after second term, 
Thus 

= -DUA^'H9,t)V{9,t)}+ (17) 
(-D^)2 {A(2)(0,i)T/(0,t)} (18) 
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where ^4*^^^ is the (fractal) drift coefficient and A^^^ is the (fractal) diffusion 
coefficient. 

Now we consider the special case with Gaussian transition probability 
P(M + .|0',<) = ^e.p{-m^i<!2)!( 



Then from equation ([121 



M„(0',i,T) = ^ / {j{e)-j{e')r 
(JW-j(r))') 

T 

The conjugate equation for moments gives 

M^{y',t,T) = ^ / (y-y')"eM- ^' }dy 

for ^^(ao) << y' and 5^(6o) >> y' we may replace the limits of the above 
integral by — oo and +oo respectively to get Mi = and M2 = t/2 For n = 1 
and n = 2 using equation (fTS]) we see that A'-^^ = 0,^^^^ = j hence the first 
term on the RHS of equation p7| vanishes and A^^'' is a constant. We can then 
write equation ([T7| as 

^V{9,t)^AiD"^)'V{9,t) (19) 

This is a new diffusion equation on the fractal curve F with "fractal" diffusion 
coefficient A. 

The ([T9| is conjugate to 

dV'(v t) 52 

Given the initial condition V'{y, 0) = S{y), the solution of the above equation 

is : 

y'(y,t) = -i=exp(-^) 
V2T:At 2At 

applying 

which gives the probability density V a,t a, location 9 — w(7i) at time t. This 
is an exact solution of the diffusion equation ([T9l) on a fractal curve F. The 
corresponding plots for the above distribution when the curve is a von-koch 
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Figure 1: Plot of (log | logV{6, t)\) against (log |^|), for a fixed t, and a straight 
line fit for it. The slope of the line is 2.4885, which is reasonably close to 2a, 
where a = log(4)/ log(3) = 1.26 is the dimension of the von Koch curve. 
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Figure 2: Plot of log{t) vs log{< >) . The slope of the line is 0.799 which is 
reasonably close to 1/a, a = log{4)/log{3) = 1.26, being the dimension of the 
von-koch curve. 
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curve in are shown in the figure [TJ from which the implicit subdiffusive 
behaviour is clear. 

In figure O the subdifFusive behaviour of motion on a fractal curve F is 
shown. F is the von-koch curve with dimension log (A) /log (3) = 1.26. The 
relation between Euclidean distance L{6) and time 't' is given by 

Jc(a,b) 

where the exponent ji decides the nature of diffusion. We find that in the above 
calculation, /x ~ 1/a, where a = 1.26, more appropriately 

^ i0.802 _ (20) 

and hence /i < 1 indicates the subdiffusive behaviour as a result of underlying 
fractal nature of space on which the particle moves. 

We conclude that the underlying fractal nature of space gives rise to subd- 
ifFusive behaviour of the diffusing entity. It is rather a deviation from gaussian 
distribution, which would have been exactly Gaussian, had the underlying space 
been ordinary and not fractal in nature. Also we see that the dimension of the 
fractal curve can be estimated from the plots of the distribution function. 

Summarizing: We have proposed a Fokker-Planck equation on fractal curves. 
An exact solution of the diffusion equation on such curves is seen to have a 
subdiffusive character. The dimension of the curve can be estimated from the 
diffusion function. 
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